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Abstract 

We present an exact solution to the non-linear equation which describes a global monopole 
in the flat space. We re-examine the metric and the geodesies outside the global monopole. 
We will see that a global monopole produces a repulsive gravitational field outside the core in 
addition to a solid angular deficit. The lensing property of the global monopole and the global 
monopole-antimonopole annihilation mechanism are studied. 

(This paper published in Class. Quantum Garv. 8 (1991) 761-767) 

1 Introduction 

The idea that magnetic monopoles ought to exist has proved to be remarkably durable. A persuasive 
argument was first put forward by Dirac [1]. Many years later, another very good argument emerged. 
The 't Hooft-Polyakov monopole solution is a synthesis of the Dirac monopole in the non-Abelian 
gauge theories [2], [3]. The monopole solutions of the realistic grand unified models based on the 
gauge groups ST/ (5) and SO(10) have been considered in [4] and [5]. An extensive list of the Kaluza- 
Klein monopole literatures was made in [6] — [11]. Surprising and qualitative phenomena arise when 
one considers the quantum mechanics of electrically charged fermions interacting with magnetic 
monopoles [12] — [16]. 

Recently, Barriola and Vilcnkin [17] have shown an approximate solution of the Einstein equa- 
tions for the metric outside a global monopole, resulting from a global symmetry breaking. Such a 
monopole has Goldstone fields with energy density decreasing with the distance only as r~ 2 , so that 
the total energy is linearly divergent at large distance. Neglecting the mass term, the monopole 
metric describes a space with a deficit solid angle. The area of a sphere of radius r is not 47rr 2 , but 
47r(l — 87rGf? 2 )r 2 . Requiring that the mass density in such a monopole should not greatly exceed the 
critical density implies there is at most one global monopole in the local group of galaxies. Equally 
stringent bounds are also derived which do not depend on cosmological assumptions, using the large 
tidal gravitational forces associated with the global monopole [18]. 

In this paper, we present an exact solution to the non-linear equation which describes the global 
monopole in the flat space. We re-examine the metric outside a global monopole. We will see a 
repulsive gravitational field outside the core in addition to a solid angular deficit. We examine the 
geodesies and find that the deflected angle is small. For ultra-relativistic particles, the deflected 
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angle also depends upon the impact parameter. If the impact parameter is of a galaxy scale and 
8 is of the grand unification scale, the term with impact parameter can be neglected. However, if 
the term is of the scale of a mini soliton star [19], this term is important. The lensing property 
of the global monopole, as well as other classical effects, are studied. Wc show that the repulsive 
gravitation force between global monopole and antimonopole (MM) does not change the MM 
annihilation proposed by Barriola and Vilenkin [17]. 



2 An exact solution of the global monopole in flat space 

The simplest model that gives rise to the global monopole is described by the Lagrangian 

jsf - \d^ a d»r - \\(r<p a - v 2 ) 2 (i) 

where (j> a is a triplet of scalar fields, a — 1,2,3. The model has a global 0(3) symmetry which is 
spontaneously broken to U(l). The field configuration describing a monopole is 

r=r]f{r) ^_ (2) 
The field equations for <fi a in the flat metric reduce to a single equation for f(r) 

■f" + -J'-^f- I{P p l) =° ( 3 ) 
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where 8 = (r^VX) -1 is the core radius of the monopole. The function f(r) grows linearly when r < 8 
and exponentially approaches unity as soon as r > 8. Barriola and Vilenkin [17] took / = 1 outside 
the core which is an approximation to the exact solution. 

On the other hand, Lan and Wang [20] have discussed an exact soliton solution of the ice-like 
structure in the form of a series of hyperbolic functions. In similar ways, we will give an exact solution 
of the global monopole. Furthermore, we also discuss the convergence of these series solutions. 

Setting x = r/5 in equation ([3]), we obtain 

x 2 f" + 2xf - x 2 f + (x 2 ~ 2)/ = (4) 
which satisfies the boundary conditions 

/(0) = , (5) 

/(oo) = 1. (6) 

This problem can be solved explicitly by the method of hyperbolic functions, in the form 

oo 

/(x) = ^c n tanh 2 ™ +1 ^ (7) 
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which satisfies the boundary condition at x — ► oo 

oo 



n=0 

By using the formula 



= V2tanh^ Vtanh 2n ^(2n + 1)- 1 , (9) 
v2 „ =0 V2 



we have 



x 2 = 2tanh 2 ^V — V-i— tanh 2 ™^ . (10) 

For the determination of c n , we first substitute equations ([7]), ([8]) and (fT0|) into equation ((4]). We 
show that the recursion formula of the coefficients c„ can be expressed as follows: 

n(2n + 3)c n = ^ ^_^+iL_ + Ak{k + 1)dl _ ( 2 fc + l)^_ x 

-(fc+l)(2fc + l)di_ijc fe + £ diCiCjCk, (11) 

' i+j+k+l=n-2 



where 

if / = -1 



(12) 



tttE^tt if / = 0,1,2,3,.... 

1=0 

Equations ([8]) and (fTTj) lead to a simple method to determine c ra . 

By using equation (|8|), we have |/(ir)| < 1 for x > 0. The series of hyperbolic functions ([7|) 
converges uniformly in the region x > 0. Thus we believe that the function series ([7]) shows an exact 
solution of equation (01 . 

Next, we consider the approximate solution of Nth order 



N 



n=0 

which satisfies the infinite boundary condition 

N 



/W=^ c Wtanh 2 "+ 1 ^= (13) 



n=0 



The approximate solution of zero order is 



/(°)=tanh^. (15) 
V2 
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The approximate solution of 10th order is 

/ 10 (a;) = 0.7539 tanh^= +0.1005 tanh 3 + 0.0414 tanh 5 -^= + 
v2 v2 v2 

+ 0.0249 tanh 7 4= + 0.0178 tanh 9 + 0.0141 tanh 11 + 
\/2 V2 \/2 

+ 0.0118 tanh 13 4= + 0.0103 tanh 15 -^= + 0.0092 tanh 17 4= + 
V2 V2 V2 

+ 0.0084 tanh 19 + 0.0078 tanh 21 (16) 
v 2 v 2 

For AT > 10, the relevant error is less than 0.5%. When N increases, the relevant error will decrease. 
One can in principle reach an arbitrary accuracy if one preserves sufficient terms in equations (|13[) 
and HU . 

In the region x ^> 1, the solution ([7|) can be written as 

f(x) = 1 - bc-^ x + o(e-^ x ) (17) 

where 

oo N 

b=Y J 2 ( 2n + 1 )c n = lim ^2(2n+l) C W. (18) 

3 The metric around a global monopole 

The most general static metric with spherical symmetry can be written as 

ds 2 = B(r)dt 2 - A(r)dr 2 - r 2 (d(9 2 + sin(9 2 d^ 2 ). (19) 

The non-vanishing components of the Ricci tensor for this metric are 

B" B' (A! B'\ IB' 

B" B' ( A' B'\ I A' 
R -=W-W{A + B)--rA> (21) 

^- i+ 2^H + l) + i' ^ 

Rtptp — sin 8 Reg. (23) 



The energy-momentum tensor of the monopole is given by 



/ ^Q+^^-l* «- ,4! 



j /2 , .. 2 /_! + i 

n 2 4 



T r r = ~r ] 2J —+r 1 2 ^ + -\r 1 \f 2 -iy , (25) 



T s e =T%= V 2 Q + h V 4 (f 2 -l) 2 . (26) 
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In the flat space the monopole core has size 8 ~ (?y\/A) _1 . For 77 <C m p , where m p is the planck 
mass, we expect that gravity does not substantially change the structure of the monopole at small 
distance so that the flat space estimation of 5 still applies. 

The useful Einstein equations are given by 

SttGT 1 t , (27) 

■ 8%GT r r . (28) 
From equation (p?Tj) , we get the general relation for A(r), 

j4 - 1 (r) = l-— /V t r 2 dr (29) 
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and from equation (|28|) . we get the general relation for B(r), 



B{r) = A _1 (r)exp(-8/i / / 2 rdr). (30) 



Substituting the solution f(r) into and (13TJ)) . we get the solutions for A(r) and -B(r). In the 
linear approximation, i.e, assuming that both A(r) and B(r) are very close to unity, the result is 

A(r) = l + 8 M +— , (31) 
r 

Cm 

B(r) = 1-8(1+^-—. (32) 

r 

where /1 = ttGij 2 , and 



m = 8^ 2 jf If 2 + ^ + ^ 2 (/ 2 - I) 2 J r 2 dr (33) 

= / f' 2 dr. (34) 
Jo 

By using equation (|16J) . we obtain the numerical result for m, 

m = -46.99?/ 2 (5. (35) 

In flat space, the core of a global monopole has size 5 ~ A - ^?? -1 . We expect that gravity does 
not notably change the structure of the monopole at small distance, so that, 

\m\ ~ 10? ? 2 <5. (36) 

There is a tiny repulsive gravitational potential due to the mass term. A freely moving particle near 
the core experiences an outward proper acceleration: 

¥= _Gm = GH (3?) 



■5 



Since |m| ~ 10r) 2 5, the mass term ~ 10rj 2 S/r. When r ^> S, the mass term may be negligible. 
Neglecting the mass term and rescaling the variables r and t, we can rewrite the metric as 

ds 2 = dt 2 - dr 2 - (1 - 8fi)r 2 (dd 2 + siTL9 2 dip 2 ). (38) 

The metric (|3"B"|) describes a space with a solid deficit angel. 



4 The geodesies 



Let us now write down the equations for the geodesies in the metric (fll?)) . From 

dV da" dx A _ 

dp 2 " A dp dp 

we have 



J (hp 

dp 



J 



dt 



dp B(r) 



(39) 

(40) 
(41) 
(42) 



Here J and E are integral constants. J represents the angular momentum of the trajectory and is 
the ratio between the proper time along the trajectory and the affine parameter p, i.e. ds 2 = — Edp 2 . 
From equations (|40|) — (|42|l thc shape of the path, r = r(ip), can be found. It is given by 



dip 



dr 



= A(r) = ^A 1 ' 2 ^) 



JB{r) 



E 
J 2 



-1/2 



(43) 



And we can parametrize the trajectory in terms of the distance of closest approach to the core, ro, 
instead of the angular momentum J. Form ^\ ro = 0, we have 



J = r (l/B(r )-E) 1/2 . 



Then equation (|4*3"| reduce to 



|g| S A(r) = ^V»(, 



B{r ) 1 - EB(r) 
B(r) l-EB(r ) 



-1/2 



(44) 



(45) 



Consider a trajectory starting from a distance L 3> ro, and with velocity V such that V 2 ^> ^t. 
Most of the deflection imprinted upon the trajectory by the gravitational filed will occur in the 
region r m ro, where 



A(r) 



ro 



1 



r [l-(r /r) 2 ] 1/2 



Gm Gm 1 
1 + 4^ + — 



1 



2r AV 2 r Q l + r /r 



(46) 
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to the first order of /x and [ijV 2 ■ Also E = 1 — V 2 to this order . The deflection angle is 



£ = 7T — Alf 



(47) 



where 



Equations gHD - (08) yield 



Ay 



'jf 



A(r)dr. 



. 1 , Gm 



(48) 



(49) 



In our order of approximation, ro coincides with the impact parameter of the trajectory. The 
analogous result to equation (j4"5|) for a gauge string is e = —Att^jl with \i the string energy per unit 
length. The monopole case differs from that of a gauge string by the dependence of e upon both 
the impact parameter and the velocity of the trajectory. 

In the case of ultrarelativistic particles (V ~ 1), the deflection angel is 



-AlTfl 



1 



Gm S 
8/j5 rg 



(50) 



where Gm/8fi5 — 0(1). We argue that light rays are deflected by various angels, dependent upon 
the impact parameter ro. If ro is a galaxy scale and S is of the order of grand unification scale, 
the second term in equation (|50|) can be neglected. In this case, the deflected angle is independent 
of the impact parameter. However, if the impact parameter is of the scale of a mini soliton star, 
the second term in equation ([50]) cannot be neglected. Such an object can be formed if a global 
monopole is swallowed by a mini soliton star. 

There is a threshold velocity Vq above which the monopole acts as a convergent lens and below 
which the monopole acts as a divergent lens. By setting e — in equation (|4T)|) one obtains 



V 2 



2(1 



AlTflS To 

Gm T 



(51) 



If ro and S are of the scale of mini soliton star, Vb is about half the speed of light. When a particle 
moves at speed Vq there is no net deflection. 



5 discussion 

The crucial question now is what is the expected density of global monopoles. Barriola and Vilenkin 
[17] proposed global monopole- antimonopolc annihilation as a possible extremely efficient mechanism 
for reducing their number. The energy of the pair (MM) is E ~ rj 2 R, where R is the MM distance. 
The attractive force acting on Mand M is F and is independent upon the distance. The repulsive 
gravitational force F lcp has maximum effect at distance R ~ S, where it is 

F rep ~ ^ ~ 10W (52) 
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Because /x ~ 10~ 6 at the typical grand unification scale 10 16 GeV, the repulsive force F rop can 
be neglected according to the large attractive force F. The large attractive force between global 
monopole and antimonopole suggest that MM annihilation is very effective and that the monopolc 
overproduction problem may not exist. 
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